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PLASMA EXPANSION IN THE GEOMETRY OF A 
COLLAPSING STAR 



R. RUFFINI, L. VITAGLIANO AND S.-S. XUE 

We describe the evolution of an electron-positron-photon plasma created by 
Sauter— Heisenberg— Euler-Schwinger mechanism around a collapsing charged star 
core in the Reissner-Nordstrom geometry external to the core, in view of the ap- 
plication in the framework of the EMBH theory for gamma ray bursts. 



1. Introduction 

In 1975, following the work on the energetics of black holes, 1 Damour and 
Ruffini 2 pointed out the existence of the vacuum polarization process a la 
Sauter-Heisenberg-Euler-Schwinger 3 ' 4 around black holes endowed with 
electromagnetic structure (EMBHs), whose electric field strength exceeds 

2 3 

the Schwinger critical value £ c = -% , where c is the speed of light, e 
and m e are electron charge and mass respectively. Damour and Ruffini 
gave reasons to believe that this process is almost reversible in the sense 
introduced by Christodoulou and Ruffini 1 and that it extracts the mass 
energy of an EMBH very efficiently: this have been proved in Ref. 5. The 
vacuum polarization process around an EMBH offered a natural mechanism 
for explaining the Gamma Ray Bursts (GRBs), just discovered at the time. 
Moreover the mechanism had a most peculiar prediction: the characteristic 
energetics of the burst should be of the order of 10 54 ergs; while nothing at 
the time was known about either the distances or the energetics of GRBs. 

More recently, after the discovery of the afterglow of GRBs and their 
cosmological distance, the idea by Damour and Ruffini has been reconsid- 
ered in Refs. 6-9 where the EMBH model for GRBs is developed. The 
evidence is now that through the observations of GRBs we are witness- 
ing the formation of an EMBH and therefore are following the process of 
gravitational collapse in real time. Even more importantly, the tremendous 
energies involved in the energetics of these sources have their origin in the 
extractable energy of black holes. 

Various models have been proposed in order to extract the rotational 
energy of black holes by processes of relativistic magnetohydrodynamics 
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(see, e.g., Ref. 10). It should be expected, however, that these processes 
are relevant over the long time scales characteristic of accretion processes. 
In the present case of GRBs a sudden mechanism appears to be at work 
on time scales of the order of few seconds or shorter and they are naturally 
explained by the vacuum polarization process introduced in Ref. 2. 

All considerations on the electric charge of stars have been tradition- 
ally directed towards the presence of a net charge on the star surface in a 
steady state condition, from the classic work by Shvartsman 11 all the way 
to the fundamental book by Punsly. 12 The charge separation can occur in 
stars endowed with rotation and magnetic field and surrounded by plasma, 
as in the case of Goldreich and Julian, 13 or in the case of absence of both 
magnetic field and rotation, the electrostatic processes can be related to the 
depth of the gravitational well, as in the treatment of Shvartsman. 11 How- 
ever, in neither case is it possible to reach the condition of the overcritical 
field needed for pair creation. 

The basic new conceptual point is that GRBs are the most violent tran- 
sient phenomenon in the universe and therefore in order to realize the con- 
dition for their occurrence, one must look at a transient phenomenon. We 
propose as a candidate the most transient phenomenon possibly occurring 
in the life of a star: the gravitational collapse. The condition for the cre- 
ation of the supercritical electromagnetic field required in the Damour and 
Ruffini work has to be achieved during the process of gravitational collapse 
which lasts less than ~ 30 seconds for a mass of 1OM and the relevant 
part of the process may be as short as 10~ 2 or even 10 -3 seconds. It is ap- 
propriate to consider a numerical example here 14 (see Fig. 1). We compare 
and contrast the gravitational collapse of a star in the two limiting cases 
in which its core of M = 3M Q and radius R = Rq is either endowed with 
rotation or with electromagnetic structure. The two possible outcomes of 
the process of gravitational collapse arc considered: cither a neutron star 
of radius of 10km or a black hole. 

In the case of rotation the core has been assumed to have a rotational 
period of ~ 15 days. For such an initial configuration we have: 

E rot ~7xlO- 12 E tot ^ \E grav \ ~6xl0~ 6 £ tot <£; bQI .~4.4xl0 5 cm. (1) 

In the collapse to a neutron star we have: 

E rot ~ 0.01E tot < \E grav \ ~ QAEtot < E bar ~ 4.4 x 10 5 cm. (2) 

The very large increase in the rotational energy is clearly due to the pro- 
cess of gravitational collapse: such a storage of rotational energy is the well 
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Figure 1. Comparing and contrasting gravitational collapse to a neutron star and to a 
black hole for a star core endowed with rotation or electromagnetic structure. Repro- 
duced from Ref. [14] with the kind permission of the author. 



known process explaining the pulsar phenomena. The collapse to a black 
hole has been estimated assuming the mass-energy formula. 1 The overall 
energetics, for the chosen set of parameters, leads to a solution correspond- 
ing to an extreme black hole, for which in principle 29% of the energy is 
extractable. 

The similar process in the electromagnetic case starts from an initial 
neutral star with a magnetosphere oppositely charged from a core with 



Let us first evaluate the amount of polarization needed in order to reach 
the above relativistic condition. Recalling that the charge to mass ratio 
of a proton is q p / (jn p ^/~G^ = 1.1 x 10 18 , it is enough to have an excess 

of one quantum of charge every 10 19 nucleons in the core of the collapsing 
star to obtain such an EMBH after the occurrence of the gravitational 
collapse. Physically this means that we are dealing with a process of charge 
segregation between the core and the outer part of the star which has the 
opposite sign of net charge in order to enforce the overall charge neutrality 
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condition. 

We then have: 

E cou ~6x KT 8 £ tot < \E grav | ~ 6 x 10- 6 E tot < E bar ~ 4.4 x 10 5 cm . (4) 

In the collapse to the neutron star configuration we have: 

E Cou ~ 0.001£ tot « |£ srat , | ~ 0.1£ tot « £ bar ~ 4.4 x 10 5 cm . (5) 

Once again, the amplification of the electromagnetic energy is due to the 
process of gravitational collapse. Again, assuming the mass-energy for- 
mula, the collapse to a black hole for the chosen set of parameters leads 
to: 

M ir = 0.9975£ tot , E Cou = 2.5 x l()- 3 E tot . (6) 

It is during such a process of gravitational collapse to an EMBH that the 
overcritical field is reached. 

The process of charge segregation between the inner core and the oppo- 
sitely charged outer shell is likely due to the combined effects of rotation 
and magnetic fields in the earliest phases of the gravitational collapse of 
the progenitor star or to a process of ionization. In the following we will 
forget about the outer shell and will treat the inner core as an electrically 
charged collapsing star. 



2. Energy extraction from a supercritical EMBH 

We know from the Christodoulou-Ruffini mass formula 1 that the mass en- 
ergy of an EMBH is the sum of the irreducible mass and the electromagnetic 
energy: 

M = M ir + J, (7) 

where Q is the charge and r + is the radius of the horizon. Moreover in 
Ref. 5 it is shown that the electromagnetic energy ^-j- is stored throughout 
the region external to the EMBH and can be extracted. If the condition 

%>£c (8) 

is fulfilled the leading extraction process is a collective process based on 
the electron-positron plasma generated by Schwinger mechanism in the 
supercritical electric field of the EMBH. 5 The condition (8) implies 



GM/c 2 



.y 1 c2-10-^<t<l (9) 
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and therefore this vacuum polarization process can occur only for an EMBH 
with mass smaller than 2 • 1O 6 M0. The electron-positron pairs are produced 
in the dyadospherc of the EMBH, 6 the spherical region whose radius r^s 
satisfies £ c = We have 



The number of particles created is 6 
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The total energy stored in the dyadosphere is 6 
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The mean energy per particle produced in the dyadosphere (E) ds 
then 
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which can be rewritten as 



r_|- 



4+ £++(!+• 



(10) 



(11) 



(12) 



E d f 
N da 



IS 



/ Ac ^ Q£ ^ 3 / Ac. A Qe fio\ 



(14) 



Such a process of vacuum polarization, occurring not at the horizon but 
in the extended dyadosphere region (r+ < r < r^) around an EMBH, has 
been observed to reach the maximum efficiency limit of 50% of the total 
mass-energy for an extreme EMBH (see e.g. Ref. 6). As discussed in Ref. 5 
the e + e~ creation process occurs at the expense of the Coulomb energy 
and docs not affect the irreducible mass, which does not depend of the 
electromagnetic energy. In this sense, 8M- 1T = and the transformation is 
fully reversible. 



3. The EMBH Theory 

In a series of papers, 6-9 Ruffini and collaborators have developed the 
EMBH theory for GRBs, which has the advantage, despite its simplicity, 
that all eras following the process of gravitational collapse to the EMBH 
are described by precise field equations which can then be numerically in- 
tegrated. 
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Starting from the vacuum polarization process a la Sauter-Heisenberg- 
Eulcr-Schwinger in the overcritical field of an EMBH first computed in 
Ref. 2, Ruffini et al. developed the dyadosphere concept. 6 

The dynamics of the e + e _ -pairs and electromagnetic radiation of the 
plasma generated in the dyadosphere propagating away from the EMBH 
in a sharp pulse (PEM pulse) has been studied by the Rome group and 
validated by the numerical codes developed at Livermore Lab. 15 ' 16 

The collision of the still optically thick e + e~-pairs and electromagnetic 
radiation plasma with the baryonic matter of the remnant of the progen- 
itor star has been again studied by the Rome group and validated by the 
Livermore Lab codes. 15,16 The further evolution of the sharp pulse of pairs, 
electromagnetic radiation and baryons (PEMB pulse) has been followed for 
increasing values of the gamma factor until the condition of transparency 
is reached. 17 

As this PEMB pulse reaches transparency the proper GRB (P-GRB) is 
emitted 8 and a pulse of accelerated baryonic matter (the ABM pulse) is in- 
jected into the interstellar medium (ISM) giving rise to an afterglow. Thus 
in GRBs we can distinguish an injector phase and a beam-target phase. 
The injector phase includes the process of gravitational collapse, the for- 
mation of the dyadosphere, as well as the PEM pulse, the cngulfmcnt of the 
baryonic matter of the remnant and the PEMB pulse. The injector phase 
terminates with the P-GRB emission. The beam-target phase addresses 
the interaction of the ABM pulse, namely the beam generated during the 
injection phase, with the ISM as the target. It gives rise to the E-APE (Ex- 
tended Afterglow Peak Emission) and the decaying part of the afterglow. 
The existence of both the P-GRB and the E-Ape is shown in Fig. 2, where 
the fit of observational data relative to GRB 991216 within the EMBH 
theory is reported. 

4. Gravitational Collapse of an Electrically Charged Core: 
Formation of Dyadosphere 

We now turn to the details of the formation of dyadosphere. If the electric 
field of a charged star core is stable against vacuum polarization during the 
gravitational collapse, 18 then an enormous amount of pairs can be created 
by Schwinger mechanism. Moreover the pairs thermalize to a positrons- 
electrons-photons plasma configuration (see Refs. 6,19, 20). Such a plasma 
undergoes a relativistic expansion. The evolution of the system and the 
details of GRB emission, along the lines summarized in the previous section, 
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Figure 2. The overall description of the EMBH theory applied to GRB 991216. The 
BATSE noise threshold is represented and the observations both of the P-GRB and of the 
E-APE arc clearly shown in the subpanels. The continuos line in the picture represents 
the theoretical prediction of the EMBH model. 



were described in Refs. 6-8, 15-17. In the latter papers the time scale of 
the gravitational collapse is neglected with respect to the hydrodynamic 
time scale. In this paper we relax this approximation: our main aim is to 
describe how the plasma expansion occurs during the gravitational collapse. 
In a forthcoming paper 18 we will discuss how the expansion is affected by 
the strong gravitational field near the horizon of the forming EMBH. 

In Refs. 20 and 21 it was suggested that the exact solution of Einstein- 
Maxwell equations describing the gravitational collapse of a thin charged 
shell can be used as an analytical model for the gravitational collapse of a 
charged core. First we briefly review some of the results of Ref. 21. The 
region of space-time external to the collasping core is Reissner-Nordstrom 
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with line element, in Schwarzschild like coordinates (t, r, 6, <ft) , 

ds 2 = -a 2 dt 2 + oT 2 dr 2 + rW, (15) 

where dtt 2 = d9 2 + sin 2 8d(f> 2 , a 2 = a 2 (r) = 1 - + M is the total 
energy of the core as measured at infinity and Q is its total charge. Let us 
label with r$ and to the radial and time- like coordinates of the shell, then the 
electromagnetic field strength on the surface of the core is £ — £ (r ) = % 
and the equation of core's collapse is 21 

^ = -^V^ 2 M-« 2 M (16) 

where H (r ) = — ^Mo'ro anc ^ ^ s rest mass °^ tne shell. The 
analytical solution of Eq. (16) was found in Ref. 21 in the form 

to = to (r ) • (17) 

Dyadosphere is formed since the instant ids = to (rd s ) when £ = E c . In the 
following we put ids = 0. 

5. Formation of e+e - Pairs around a Collapsing Charged 
Core 

For t < ids, £ < £ c and the Schwingcr process of e+e~ pairs creation is ex- 
ponentially suppressed. For t > ids the Schwinger process becomes relevant 
and e + e~ pairs are created. As shown in Refs. 19,20 the pairs created at 
radius ro < rd s oscillate with ultrarelativistic velocity and partially annihi- 
late into photons. At the same time the electric field oscillates around zero 
and the amplitude of such oscillations decreases in time: in a time of the 
order of 10 2 — 10 4 H/m e c 2 the electric field is effectively screened to about 
the critical value; more precisely, the average of the electric field £ over one 
period of oscillation is 0, but the average of £ 2 is of the order of £ 2 . As a 
result an energy density has been deposited 5 on the pairs and the photons 
given by 



1 



eo(ro) = -[£ 2 (r ) -£?] = ,„ 



(18) 



The pairs and the photons are expected to thermalize 6 ' 19 ' 20 , to an e+e 7 
plasma equilibrium configuration: 

n e + = n e - ~ n 7 = n , (19) 
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(where n, is the proper number density of particles of type •), and reach 
an average temperature To such that 

e (T ) ee e 7 (T ) + e e+ (T ) + e e - (T ) = e ; (20) 

here e. (T) is the equilibrium proper energy density at temperature T for 
the species •. Then n e ± (n 7 ) are given by Fermi (Bose) integrals once the 
temperature T is known. 

6. Plasma's Expansion 

The highly energetic plasma so formed undergoes a relativistic expansion. 
As will be shown, the expansion (hydrodynamic) time-scale is much bigger 
than both the pair creation and the thermalization time-scales, then the 
process can be described as follows: at any time to it begins to expand a 
slab of plasma of thickness Al = a _1 Ar (as measured by static observers) 
produced at radius r = r (to)- Al can be chosen very small in comparison 
with rds so that, in particular, the temperature T is approximately constant 
in the slab. Moreover Al has to be much bigger than the quantum length 
scale (~ h/m e c). 

We can follow the expansion of each slab of plasma by using conservation 
of energy and number of particles. Note that Eqs. (20) and (19) provide 
initial data for the problem of the expansion. We describe the expansion 
of a single slab using the following approximations: 

(1) the geometry in which the expansion occurs is Reissner- Nordstrom 
with line-element given by (15). In particular we will denote by £ a 
the static vector field normalized at unity at spatial infinity, and 
by {£t} t the family of space- like hypersurfaces orthogonal to £ a (t 
being the Killing time); 

(2) the plasma is assumed to be a neutral perfect fluid characterized 
by proper energy density e, proper pressure p, proper particle (elec- 
trons, positrons and photons) number density n and 4— velocity tt°; 

(3) at any instant, electrons, positrons and photons in a single slab are 
assumed to be at thermal equilibrium with temperature T, possibly 
different from slab to slab. The slabs are uncorrelated in the sense 
that they do not share energy nor particles. In other words the 
expansion of each slab is adiabatic; this will be checked a posteriori 
(see also Ref. 16); 

(4) the thickness Al = oT 1 Ar of a slab as measured by static observers 
is constant. 
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Figure 3. A sequence of snapshots of coordinate baryon energy density is shown from 
the numerical solution of partial differential continuity equations. This run correspond 
to an EMBH of mass M = 1O 3 M and charge to mass ratio £ = Q/M = 0.1. 



The last approximation is justified by the result in Ref. 16, where it is 
shown, by numerical integration of the partial differential continuity equa- 
tions, that the baryon energy density of an expanding slab of plasma en- 
riched with nucleons from the remnant of the progenitor star is localized in 
a region of constant thickness (see Fig. 3). 

Given the above assumptions, both the energy momentum-tensor T ab = 
(e + p) u a u b + pg ab and the electron (positron) -number current n° = n e u a 
are conserved: 

X7 b T ab = 0, (21) 
V b n b = 0. (22) 

In particular, using assumption (4) one can reduce the partial differential 
continuity equations (21) and (22) to ordinary differential equations for the 
radial coordinate r and the temperature T of the single slab as functions 
of time (see Ref. 16). The equation of motion of a single slab can be 
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numerically integrated with initial conditions 



r(to) 



(23) 



dr I 

dt \t=t { 



(24) 
(25) 



T(to) 



T . 



The overall motion of the plasma is the superposition of motions of single 
shells. The typical plasma expansion curves are shown in Fig. 4 from the 
numerical integration of the equations of motion. 

The curvature of space-time strongly affects the motion of plasma in 
the vicinity of the EMBH horizon and in turn the phenomenology of the 
GRB. We discuss these issues in a forthcoming paper. 18 
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Figure 4. Collapse curve of an elettromagnetic star core (dashed line) as derived by 
Eq. (16) in the case M = Mo = 20Mq, Q = 0.1M and expansion curves of plasma. 



